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FOREWORD 

The  range  tables  for  Navy  guns  include  a  nomogram  for  computing  ballistic 
density,  which  may  be  used  when  only  the  surface  density  is  known.  While  these 
nomograms  are  intended  for  use  when  better  estimates  of  ballistic  density  are  not 
available,  they  are  based  on  data  prepared  in  1922,  despite  the  advent  of  radio¬ 
sondes  and  the  establishment  of  ocean  weather  stations. 

Since  there  frequently  are  times  when  a  radiosonde  observation  cannot  be  made 
in  the  immediate  operating  area  just  befojre  firing,  or  a  timely  and  applicable 
meteorological  message  from  the  Fleet  Numerical  Weather  Center  (FNWC)  is  not  avail¬ 
able,  another  method  for  obtaining  ballistic  densities  would  serve  a  useful 
purpose . 


This  report  documents  a  study  that  shows  that  ballistic  densities  can  be 
estimated  from  the  surface  density  with  a  reasonable  degree  of  accuracy.  Data  are 
also  presented  about  the  upper  atmosphere  which  may  be  useful  when  further  studies 
pertaining  to  the  atmosphere  are  conducted. 

The  study  was  funded  by  the  Naval  Sea  Systems  Command  Surface  Weapons 
Aerodynamics  and  Structures  Research  and  Development  Block  Program. 


Released  by: 

0.  F.  BRAXTON,  Head 
Strategic  Systems  Department 
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INTRODUCTION 

One  of  the  basic  premises  of  fire  control  systems  in  use  in  the  Fleet  today  is 
a  "Standard  Atmosphere."  Obviously,  the  same  premise  is  also  used  in  the  prepara¬ 
tion  of  ballistic  tables  for  the  Fleet.  Since  the  conditions  used  as  a  basis  for  a 
standard  atmosphere  are  rarely  encountered  on  earth,  corrections  for  nonstandard 
conditions  are  always  in  order.  A  procedure  frequently  used  to  correct  for  a  non¬ 
standard  density  involves  a  ballistic  air  density  in  the  solution  of  fire  control 
problems.  The  ballistic  air  density  may  be  defined  as  a  scale  factor  that,  when 
multiplied  by  the  standard  air  density  at  each  altitude,  produces  a  density  struc¬ 
ture  that  would  cause  the  same  displacement  of  a  projectile  impact  point  as  the 
actual  density  structure. 

The  most  reliable  method  of  determining  ballistic  density  is  to  use  accurate, 
representative,  and  timely  meteorological  observations  in  the  computation.  There 
are,  however,  occasions  when  these  measurements  are  not  available,  and  alternate 
methods  are  necessary  to  provide  the  requisite  information.  It  is  a  standard  pro¬ 
cedure  to  include  a  nomogram  in  the  ballistic  tables  for  the  Fleet.  The  use  of  the 
nomogram  is  an  alternative  that  is  available  to  the  Fleet  for  estimating  the  upper 
air  density  for  naval  gunfire  when  a  ballistic  meteorological  forecast  is  not 
available.  The  nomograms  may  be  used  to  estimate  the  ballistic  density  when  the 
air  temperature,  pressure,  and  relative  humidity  at  the  surface  are  known. 

The  basic  data  used  to  prepare  the  nomograms  were  obtained  from  Table  IV,  The 
Ballistic  Density  Factor,*  which  was  published  in  1935.  A  footnote  to  this  table 
states,  "This  table  has  been  extracted  from  Table  C  of  a  pamphlet  entitled  Method 
for  Determining  Air  Temperature  and  Ballistic  Air  Density,  published  by  the  Signal 
Corps  (Meteorological  Service),  U.S.  Army  in  1922." 

Fire  control  systems  have  been  improved  to  the  point  diere  an  accurate  deter¬ 
mination  of  the  upper  air  density  is  required  if  the  potential  of  the  systems  is  to 
be  fully  utilized.  The  need  for  improving  the  nomograms  or  providing  an  alternate 
procedure  to  determine  the  ballistic  density  has  long  been  recognized. 

The  work  described  in  this  report,  stated  briefly,  consisted  of  (1)  obtaining 
Rawinsonde  data  from  various  weather  stations,  (2)  sorting  the  data  according  to 
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surface  densities  and  seasons,  (3)  fitting  the  upper  air  densities  with  a  least- 
squares  technique,  and  (4)  computing  ballistic  densities. 

Rawinsonde  data  used  in  the  study  were  obtained  from  five  maritime  stations 
which  took  soundings  daily  at  0000  and  1200  hr  G.m.t.  over  a  period  from  1968  to 
1973.  The  data,  consisting  of  both  significant  and  mandatory  levels,  were  provided 
by  the  National  Climatic  Center,  Asheville,  North  Carolina.  Mandatory  levels  are 
specified  pressure  altitudes,  which  are  obtained  by  interpolation.  Significant 
levels  are  the  altitudes  that  are  considered  essential  in  defining  a  weather  pro¬ 
file.  Soundings  that  did  not  include  surface  density  were  not  used  in  the  study. 
Some  of  the  obvious  errors  were  eliminated  in  processing  the  data  prior  to  perform¬ 
ing  the  least-squares  fit,  but  no  attempt  was  made  to  eliminate  errors  in  the  basic 
data  beyond  tt:ai.  point. 

The  geodetic  coordinates  of  the  five  stations  used  in  the  study  are  listed 
below: 


Station  Name 

Latitude 

Longitude 

E 

35°  O'  N 

48°  0*  W 

V 

34°  O’  N 

164’  O’  E 

C 

52°  45’  N 

35’  30*  W 

D 

44’  O’  N 

41’  O’  W 

N 

30°  O’  N 

140’  0’W 

A  map  of  the  stations  is  shown  in  Figure  1. 

No  land-based  stations  were  used  in  the  analysis,  even  though  naval  gunfire 
support  is  conducted  on  beach  areas.  Since  the  density  correlations  are  markedly 
affected  by  the  nature  of  the  terrain, 2  there  is  more  than  a  probability  that 
somewhat  different  results  would  be  obtained  for  different  littoral  areas  of  the 
earth. 

All  computer  programs  used  in  the  study  were  coded  in  FORTRAN  EXTENDED  and 
were  executed  with  a  Control  Data  Corporation  (CDC)  6700  computer. 

The  purpose  of  performing  the  work  described  in  this  report  was  to  establish  a 
procedure  to  determine  the  correlation  between  upper  air  and  surface  densities  and 
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to  devise  a  method  for  estimating  the  ballistic  air  density  based  on  surface 
conditions . 


DETAILED  PROCEDURE 


The  first  step  in  the  procedure  was  to  compute  the  air  density  using  the  tem¬ 
perature,  pressure,  and  relative  humidity.  The  geopotential  altitudes  were  con¬ 
verted  to  geometric  altitudes.  Sondes  that  did  not  include  surface  conditions  were 
omitted.  The  data  for  each  station  were  sorted  into  two  groups  according  to  the 
time  of  day  and  stored  in  permanent  files  with  the  Indexed-Sequential  File 
Organization.  This  procedure  permits  accessing  the  data  in  either  a  sequential  or 
random  order.  The  sondes  were  assigned  indices  which  were  stored  in  separate 
sequential  files.  One  file  of  indices  was  created  with  the  corresponding  surface 
densities  in  ascending  numerical  order,  and  a  second  file  was  created  with  the  sur¬ 
face  densities  also  sorted  into  the  four  seasons.  A  follow-up  program  was  designed 
to  skip  specified  segments  when  the  files  were  read.  This  made  it  convenient  to 
group  the  data  in  various  ways  without  rearranging  the  Indexed-Sequential  files. 

The  data  were  then  selected  in  groups  where  the  median  surface  density  of  each 
group  varied  by  a  percentage  point  from  the  U.S.  Standard  surface  density  as 
defined  in  Reference  3  (i.e.,  into  classes  by  surface  density  ratios  vrtiere  the 
class  interval  was  one  percent).  Since  large  variations  were  encountered  in  the 
frequency  of  the  surface  densities,  it  was  necessary  to  eliminate  some  of  the 
sondes  in  the  mid-ranges.  Several  fits  were  made  using  various  numbers  of  samples 
in  each  fit.  No  significant  differences  were  obtained  by  using  more  than  25  sam¬ 
ples;  consequently,  a  further  restriction  was  placed  on  each  group  that  limited  the 
number  of  sondes  to  25.  While  some  of  the  groups  outside  the  mid-ranges  had  varia¬ 
tions  in  surface  densities  of  one  percent,  the  actual  number  of  samples  was,  in 
some  cases,  less  than  25.  In  contrast,  some  of  the  groups  of  25  samples  in  the 
mid-ranges  had  variations  in  surface  densities  of  only  a  fraction  of  one  percent. 

Various  functions  were  investigated  to  represent  a  mean  density  profile  for 
each  group  of  data,  Comparisons  were  made  of  three  different  functions  for  rep¬ 
resenting  the  density  in  the  altitude  range  from  0  to  36,089  ft: 
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1.  a  third-degree  polynomial 

2.  a  third-degree  polynomial  for  the  ratio  of  the  observed  density  to  the 
U.S.  Standard  density 

3.  an  exponential  function  plus  a  linear  term 

While  these  three  functions  produced  similar  results,  the  exponential  function  was 
selected.  The  fitting  technique  is  described  in  Appendix  A.  The  following  expres¬ 
sion  was  used  in  the  altitude  range  from  0  to  36,089  ft: 

Cc  h) 

p  =  cje  +  c3h  (1) 


wnere , 

p  is  the  density  (lb/ft^) 
h  is  the  altitude  (ft) 

c^,  c^,  and  c^  are  constants  determined  by  the  fit 

A  similar  expression  was  used  to  fit  the  data  in  the  region  from  36,089-  to 
65,000-ft  altitude: 


P 


(c  *(h  -  36089)) 

P  e  4 

36089 


+  c5*(h  -  36089) 


(2) 


where , 


P  36089  density  at  36,089  ft,  computed  using  equation  (1)  and  the 

data  up  to  36,089  ft 

c,  and  c,  are  coefficients  determined  in  the  fit  above  36,089  ft 
4  5 


A  bias  was  noted  in  the  lower  altitudes  with  each  of  the  three  methods.  After 
some  experimenting,  data  below  3,500  ft  were  excluded  from  the  fit  to  determine  the 
coefficients.  The  fit  was  then  extrapolated  to  mean  sea  level,  and  a  linear  cor¬ 
rection  was  added  from  0  to  3,500  ft,  which  made  P  equal  to  Pq  at  the  surface: 
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(c  h) 

P  =  +  c3h  +  (PQ  -  0^(3500  -  h)/3500  (3) 


where , 

P  is  the  density  (lb/ft^) 

Cj ,  c^,  and  are  the  coefficients  obtained  with  equation  (1) 

h  is  the  altitude  (ft) 

p0  is  the  median  of  the  observed  surface  densities  (lb/ft^) 

Although  this  procedure  does  not  use  any  of  the  available  data  between  surface 
and  3,500  ft  in  the  least-squares  fit  to  determine  the  coefficients  c^  ,  c^,  and 
c^  ,  residuals  from  the  fit  in  this  region  are  included  in  the  analysis.  Adding  the 
linear  correction  made  a  significant  improvement  in  the  residuals  from  surface  to 
3,500  ft. 


The  upper  boundary  of  36,089  ft  was  selected  from  the  zone  where  equation  (1) 
was  used  in  order  to  coincide  with  the  breakpoint  in  the  U.S.  Standard  atmosphere. 
While  there  are  some  seasonal,  latitudinal,  and  diurnal  variations  in  the  isopycnic 
level  where  the  density  remains  nearly  constant,  no  attempt  was  made  to  explore 
this  facet . 

After  fitting  the  observed  density  data,  using  equations  (1)  and  (2),  density 
profiles  were  computed  using  equations  (2)  and  (3)  with  the  coefficients  obtained 
in  the  least-squares  fit.  Ballistic  densities  for  each  profile  were  computed  with 
the  equation. 


n 

Ballistic  Density  =  £  (WF.  *  R.) 

i=  1  1  1 


(4) 


where , 


n  is  the  number  of  zones 

WFj  is  the  density  weighting  factor  for  the  ith  zone 
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is  the  ratio  of  the  mean  density  of  the  ith  zone  of  the  profile  to  the  U.S. 

Standard  density  at  that  altitude,  which  is  the  mid-point  of  the  ith  zone. 

The  density  weighting  factors  used  to  compute  the  ballistic  density  were 
obtained  by  computing  trajectories  with  ballistic  parameters  for  the  5-Inch  High 
Fragmentation  Projectile  fired  with  a  Mk  73  charge,  which  produces  an  initial 
velocity  of  2,950  ft/ sec.  The  altitudes  covered  by  each  zone  are  presented  in 
Table  1,  and  a  tabulation  of  the  weighting  factors  are  presented  in  Table  2. 


TABLE  1.  ALTITUDE  OF  ZONES 


Zone 

Altitude  (ft) 

Zone 

Altitude  (ft) 

Zone 

Altitude  (ft) 

1 

0 

_ 

1,000 

10 

9,000 

_ 

10,000 

19 

26,000 

28,000 

2 

1,000 

- 

2,000 

11 

10,000 

- 

12,000 

20 

28,000 

- 

30,000 

3 

2,000 

- 

3,000 

12 

12,000 

- 

14,000 

21 

30,000 

- 

32,000 

4 

3,000 

- 

4,000 

13 

14,000 

- 

16,000 

22 

32,000 

- 

34,000 

5 

4,000 

- 

5,000 

14 

16,000 

_ 

18,000 

23 

34,000 

— 

36,000 

6 

5,000 

- 

6,000 

]  5 

18,000 

- 

20,000 

24 

36,000 

- 

38,000 

7 

6,000 

- 

7,000 

16 

20,000 

- 

22,000 

25 

38,000 

- 

40,000 

8 

7,000 

- 

8,000 

17 

22,000 

- 

24,000 

26 

40,000 

- 

45,000 

9 

8,000 

- 

9,000 

18 

24,000 

- 

26,000  , 

27 

45,000 

- 

50,000 

FITTING  A  SAMPLE  SET  OF  DATA  WITH  AN  ALTERNATE  PROCEDURE 

An  alternate  method  was  used  to  fit  a  sample  set  of  data  in  the  region  from 
surface  to  36,000  ft.  In  this  procedure,  the  virtual  temperature  was  computed, 
using  the  observed  temperature  and  relative  humidity.  Separate  least-squares  fits 
were  then  made  to  the  virtual  temperatures  and  pressures.  The  temperature  and 
pressure  profiles  obtained  in  this  manner  were  then  used  to  compute  the  density. 
This  procedure  and  the  results  are  described  in  detail  in  Appendix  B.  The  sample 
set  of  data  was  obtained  from  25  sondes  that  had  been  selected  previously  for  fit¬ 
ting  with  equation  (1).  The  sondes  were  observed  at  Station  E  during  the  spring 
season,  1200  hr  G.m.t.,  on  days  that  had  a  surface  density  close  to  the  U.S. 
Standard  density  at  mean  sea  level. 
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The  purpose  of  using  Che  alternate  method  was  to  evaluate  the  method  described 
in  the  section  Procedure  and  to  provide  additional  information  about  the  atmo¬ 
sphere,  particularly  in  the  lower  altitudes. 

While  this  method  is  a  good  procedure  for  providing  temperature,  pressure,  and 
density  profiles,  it  is  a  more  time-consuming  method  for  obtaining  density  profiles 
than  the  method  of  making  a  least-squares  fit  to  the  densities. 

The  density  profile  obtained  with  this  method  did  not  result  in  a  noticeably 
different  profile,  but  it  did  provide  additional  information  about  the  temperature 
and  pressure.  The  pressure  data  from  surface  to  36,000-ft  altitude  were  fit  with  a 
single  function  that  has  a  continuous  derivative.  The  temperature  data  were 
divided  into  three  altitude  zones.  Separate  least-squares  fits  were  made  to  each 
zone  in  a  manner  that  made  the  curve  continuous  but  with  a  different  lapse  rate  in 
each  zone.  The  improvement  in  using  the  correction  term  in  equation  (3)  to  fit  the 
density  is  largely  due  to  the  variations  in  the  temperature  lapse  rate  at  the  lower 
altitudes . 


FITTING  THE  U.S.  STANDARD  ATMOSPHERE 


In  order  to  further  test  the  least-squares  method  and  the  exponential  repre¬ 
sentation  of  the  density,  a  fit  was  made  to  the  U.S.  Standard  density,  using  equa¬ 
tions  (1)  and  (2).  Since  the  U.S.  Standard  temperature  has  a  constant  lapse  rate 
from  surface  to  36,089-ft  altitude,  the  correction  term  in  equation  (3)  was  not 
used.  It  should  be  noted  that  equations  (1)  and  (2)  were  not  intended  to  define 
the  U.S.  Standard  density  but  were  selected  to  fit  the  observed  data.  However,  the 
standard  deviations  in  the  fit  were  only  0.14  and  0.28  percent  in  the  zones  from  0- 
to  36,089-ft,  and  36,089-  to  65,000-ft  altitude,  respectively.  Since  the  density 
varies  with  altitude,  the  residuals  were  expressed  as  ratios  of  the  U.S.  Standard 
density  at  the  corresponding  altitude.  A  tabulation  of  the  fit  is  presented  in 
Tables  3  and  4. 
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TABLE  3.  LEAST-SQUARES  FIT  OF  U.S.  STANDARD  DENSITY 

0  -  36,000  ft 


1 

2 

3 

4 

U.S.  Standard 

L.  S.  Fit 

Col.  2  -  Col.  3 
Col.  3 

Altitude 

Density 

Density 

(ft) 

(lb/ ft  ) 

( lb/ ft  ) 

0. 

0.076474 

0.07665845 

-0.00240610 

500. 

0.075362 

0.07551388 

-0.00201122 

1,000. 

0.074262 

0.07438439 

-0.00164533 

1,500. 

0.073174 

0.07326977 

-0.00130712 

2,000. 

0.072098 

0.07216982 

-0.00099517 

2,500. 

0.071035 

0.07108433 

-0.00069395 

3,000. 

0.069984 

0.07001309 

-0.00041553 

3,500. 

0.068945 

0.06895591 

-0.00015824 

4,000. 

0.067917 

0.06791259 

0.00006494 

4,500. 

0.066902 

0.06688293 

0.00028510 

5,000. 

0.065898 

0.06586675 

0.00047450 

5,500. 

0.064906 

0.06486384 

0.00064991 

6,000. 

0.063925 

0.06387404 

0.00079783 

6,500. 

0.062956 

0.06289715 

0.00093568 

7,000. 

0.061998 

0.06193299 

0.00104969 

7,500. 

0.061051 

0.06098138 

0.00114159 

8,000. 

0.060116 

0.06004216 

0.00122986 

8,500. 

0.059191 

0.05911513 

0.00128338 

9,000. 

0.058278 

0.05820014 

0.00133776 

9,500. 

0.057375 

0.05729702 

0.00136106 

10,000. 

0.056483 

0.05640559 

0.00137245 

10,500. 

0.055602 

0.05552569 

0.00137431 

11,000. 

0.054731 

0.05465717 

0.00135085 

11,500. 

0.053871 

0.05379985 

0.00132241 

12,000. 

0.053022 

0.05295360 

0.00129174 

12,500. 

0.052182 

0.05211824 

0.00122336 

13,000. 

0.051353 

0.05129363 

0.00115745 

13,500. 

0.050534 

0.05047962 

0.00107734 

14,000. 

0.049725 

0.04967605 

0.00098539 

14,500. 

0.048926 

0.04888278 

0.00088407 

15,000. 

0.048137 

0.04809968 

0.00077598 

15,500. 

0.047358 

0.04732658 

0.00066389 

16,000. 

0.046589 

0.04656336 

0.00055067 

16,500. 

0.045829 

0.04580987 

0.00041754 

17,000. 

0.045079 

0.04506598 

0.00028881 

17,500. 

0.044338 

0.04433156 

0.00014526 

18,000. 

0.043606 

0.04360647 

-0.00001072 

18,500. 

0.042884 

0.04289057 

-0.00015329 
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TABLE  3.  LEAST-SQUARES  FIT  OF  U.S.  STANDARD  DENSITY 
0  -  36.000  ft  (Continued) 


1 

Altitude 

(ft) 

2 

U.S.  Standard 
Density 

(lb/ft3) 

19,000. 

0.042171 

19,500. 

0.041468 

20,000. 

0.040773 

20,500. 

0.040087 

21,000. 

0.039410 

21,500. 

0.038742 

22,000. 

0.038083 

22,500. 

0.037432 

23,000. 

0.036790 

23,500. 

0.036156 

24,000. 

0.035531 

24,500. 

0.034914 

25,000. 

0.034306 

25,500. 

0.033705 

26,000. 

0.033113 

26,500. 

0.032529 

27,000. 

0.031952 

27,500. 

0.031384 

28,000. 

0.030823 

28,500. 

0.030270 

29,000. 

0.029725 

29,500. 

0.029187 

30,000. 

0.028657 

30,500. 

0.028134 

31,000. 

0.027619 

31,500. 

0.027110 

32,000. 

0.026610 

32,500. 

0.026116 

33,000. 

0.025629 

33,500. 

0.025149 

34,000. 

0.024676 

34,500. 

0.024210 

35,000. 

0.023751 

36,000. 

0.022853 

Mean  of  Column  4 


3 

4 

L.  S.  Fit 

Col.  2  -  Col.  3 

Density 

V 

Col.  3 

(lb/ft3) 


0.04218375 

0.04148588 

-0.00030232 

-0.00043087 

0.04079682 

0.04011645 

0.03944466 

0.03878132 

0.03812631 

-0.00058376 

-0.00073413 

-0.00087867 

-0.00101387 

-0.00113606 

0.03747953 

0.03684084 

0.03621014 

0.03558732 

0.03497226 

-0.00126803 

-0.00137995 

-0.00149516 

-0.00158251 

-0.00166590 

0.03436486 

0.03J76501 

0.03317260 

0.03258753 

0.03200970 

-0.00171279 

-0.00177724 

-0.00179668 

-0.00179614 

-0.00180254 

0.03143900 

0.03087534 

0.03031861 

0.02976872 

0.02922557 

-0.00174942 

-0.00169508 

-0.00160323 

-0.00146858 

-0.00131974 

0.02868907 

0.02815913 

0.02763564 

0.02711853 

0.02660771 

-0.00111786 

-0.00089227 

-0.00060224 

-0.00031466 

0.00008623 

0.02610307 

0.02560455 

0.02511204 

0.02462548 

0.02414477 

0.00049524 

0.00095501 

0.00147165 

0.00205161 

0.00270177 

0.02366983 

0.02273694 

0.00342940 

0.00510442 

-0.00000178 


Standard  Deviation  of  Column  4  0.00140643 
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TABLE  4.  LEAST-SQUARES  FIT  OF  U.S.  STANDARD  DENSITY 
36,200  -  65,000  ft 


1 

2 

3 

4 

U.S.  Standard 

L.  S.  Fit 

Col.  2  -  Col.  3 
Col.  3 

Altitude 

Density 

Density 

o 

(ft) 

( lb/ ft; 

( lb/ft  ; 

36200. 

0.022666 

0.02253574 

0.00578032 

37000. 

0.021814 

0.02169457 

0.00550490 

38000. 

0.020794 

0.02068713 

0.00516582 

39000. 

0.019822 

0.01972648 

0.00484233 

40000. 

0.018895 

0.01881043 

0.00449577 

41000. 

0.018012 

0.01793693 

0.0041854 1 

42000. 

0.017170 

0.01710398 

0.00385968 

43000. 

0.016367 

0.01630972 

0.00351190 

44000. 

0.015602 

0.01555234 

0.00319287 

45000. 

0.014873 

0.01483014 

0.00289033 

46000. 

0.014178 

0.01414147 

0.00258344 

47000. 

0.013516 

0.01348478 

0.00231541 

48000. 

0.012884 

0.01285858 

0.00197664 

49000. 

0.012282 

0.01226147 

0.00167447 

50000. 

0.011709 

0.01169208 

0.00144691 

51000. 

0.011162 

0.01114914 

0.00115366 

52000. 

0.010641 

0.01063141 

0.00090240 

53000. 

0.010144 

0.01013772 

0.00061976 

54000. 

0.0096701 

0.00966695 

0.00032545 

55000. 

0.0092186 

0.00921805 

0.00005944 

56000. 

0.0087882 

0.00879000 

-0.00020437 

57000. 

0.0083780 

0.00838182 

-0.00045561 

58000. 

0.0079870 

0.00799260 

-0.00070018 

59000. 

0.0076142 

0.00762145 

-0.00095105 

60000. 

0.0072589 

0.00726754 

-0.00118830 

61000. 

0.0069202 

0.00693006 

-0.00142260 

62000. 

0.0065973 

0.00660825 

-0.00165751 

63000. 

0.0062895 

0.00630139 

-0.00188719 

64000. 

0.0059961 

0.00600878 

-0.00211035 

65000. 

0.0057164 

0.00572976 

-0.00233130 

Mean  of  Column  4  0.0000 3016 

Standard  Deviation  of  Column  4  0.00284086 

The  following  equations  were  used  in  the  least-squares  fit  of  the  U.S. 
Standard  density: 
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0  <  h  i  36089 

(c2h) 

p  ■  c^e  +  c^h 


36,089  <  h  i  65,000 

(c  *  (h  -  36089) 

»  -  „3608,  .  “=s  *  <h  -  36089) 


where 


h  is  the  altitude  (ft) 


The  coefficients  obtained  in  the  least-squares  fit  are: 


C2 

c3 

p 36089 

C4 

C5 


0.0766584476  (lb/ft3) 
-0.282539693789  E-4  (ft-1) 
-0.138466611596  E-6  (lb/ ft4) 
0.0226549960  (lb/ft3) 
-0.475504250122  E-4  (ft-1) 
0.459663274698  E-ll  (lb/ft3) 


RESULTS 

The  coefficients  that  were  obtained  by  fitting  the  density  data,  and  the  root 
mean  squares  (rms)  or  standard  deviation  of  the  residuals  and  their  algebraic  means 
are  given  in  Appendix  C  for  the  five  stations  analyzed.  Data  are  presented  for  two 
times  of  day  (0000  and  1200  hr  G.m.t.)  for  the  four  seasons  and  for  the  seasons 
combined  for  each  station. 

Profiles  of  the  density  ratios  for  station  E  are  given  in  Appendix  D.  The 
ratios  were  obtained  by  evaluating  the  least-squares  fit  of  the  observed  densities 
at  specified  altitudes  and  dividing  the  density  in  each  case  by  the  U.S.  Standard 
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density  for  the  same  altitude.  Profiles  are  shown  for  the  two  times  of  day  for 
each  of  the  four  seasons  and  for  the  seasons  combined. 

Profiles  of  the  ballistic  densities  obtained  by  multiplying  the  ratios  with 
the  density  weighting  factors  using  equation  (4)  are  shown  in  Appendix  E  for  sta¬ 
tion  E  for  the  same  conditions.  Density  ratios  and  ballistic  densities  were  com¬ 
puted  for  each  of  the  five  stations,  but  the  profiles  are  shown  only  for  station  E 
in  Appendixes  D  and  E. 

While  the  ballistic  densities  were  computed  for  the  27  zones  listed  in 
Table  1,  only  an  abbreviated  tabulation  is  given  in  Appendix  F.  Data  are  presented 
for  only  two  zones,  30,000  and  50,000  ft.  The  ballistic  densities  were  computed 
with  equation  (4),  using  the  density  weighting  factors  given  in  Table  2.  Data  for 
the  two  zones  are  given  for  each  of  the  five  stations  and  for  each  surface  condi¬ 
tion  where  sufficient  data  were  available  to  make  a  least-squares  fit  to  the 
observed  data. 


DISCUSSION  OF  RESULTS 

One  question  is  of  primary  interest— is  there  sufficient  correlation  between 
surface  density  and  that  of  the  upper  atmosphere  to  accurately  predict  the  bal¬ 
listic  density  for  various  altitudes  from  surface  measurements?  To  answer  this 
question,  we  first  note  that  the  density  ratio  curves  given  in  Appendix  D  for  sta¬ 
tion  E  are  families  of  very  similiar  curves.  If  some  of  the  curves  were  missing, 
we  could  do  a  fair  job  of  reproducing  the  curves  by  interpolation.  This  would  seem 
to  indicate,  qualitatively  at  least,  that  one  could  predict  the  density  ratios 
aloft  from  the  starting  point  of  the  curve  on  the  x-axis,  which  is  the  surface 
density  ratio.  The  same  is  true  of  the  curves  given  for  ballistic  density  in 
Appendix  E.  It  deserves  mentioning  that  the  points  plotted  in  Appendix  D  are  not 
the  mean  values  of  individual  sondes  but  are  values  corresponding  to  the  functional 
fit . 


The  sondes  were  grouped  in  order  that  a  more  reliable  mean  and  standard  devia¬ 
tion  could  be  established.  It  is  recognized  that  the  curves  in  Appendix  D  do  not 
reflect  the  errors  of  the  means  or  the  variations  in  the  upper  atmosphere  within  a 
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group  of  sondes.  A  total  of  14,800  sondes  were  processed  in  this  study  in  an 
attempt  to  establish  reliable  means  and  standard  deviations.  In  every  fit,  the 
mean  and  standard  deviation  were  computed.  These  are  tabulated  in  Appendix  C.  A 
large  majority  of  the  means  are  only  a  small  fraction  of  one  percent. 

The  standard  deviations  of  the  residuals  are  given  in  Table  5  according  to  percen¬ 
tiles.  A  total  of  300  fits  was  made  in  each  of  the  two  altitude  zones,  0  to  36,000 
and  36,000  to  65,000  ft.  The  standard  deviations  of  the  residuals  were  arranged  in 
numerical  order  without  regard  to  station,  season,  or  time  of  day.  The  standard 
deviation  of  the  50th  percentile  is  1.114  percent  in  the  lower  altitude  zone. 


TABLE  5.  STANDARD  DEVIATIONS  OF  THE  RESIDUALS 


Percent  Lie 

A1 1  i  t  ude 

0  -  36,000  (ft) 

<%) 

Altitude 
36,000  -  65,000 
(%) 

10 

0.675 

2.019 

20 

0.799 

2.293 

30 

0.908 

2.586 

40 

1.000 

2.832 

50 

1.114 

3.067 

60 

1.227 

3.295 

70 

1.424 

3.584 

80 

1.577 

3.959 

90 

1.703 

4.298 

100 

2.257 

6.377 

Note:  The  standard  deviations  in  this  table  were  obtained  by  multiplying  the 
standard  deviations  (SiCMA)  in  Appendix  C  by  100  in  order  to  express  them  as  a 
percentage  instead  of  as  a  ratio. 


It  should  also  be  remembered  that  there  are  variations  in  the  surface  densi¬ 
ties  because  of  grouping  the  data,  which  probably  has  a  neglible  effect  on  the  means 
but  makes  a  small  contribution  to  the  standard  deviations.  The  noise  level  in  the 
data  appears  to  increase  with  altitude,  which  is  to  be  expected.  The  residuals  are 
much  larger  in  the  fits  above  36,000  ft  than  they  are  at  the  lower  altitudes,  while 
the  density  should  be  more  stable  at  the  higher  altitudes. 


It  may  be  noted  that  the  ballistic  density  represents  a  weighted  average  of 
the  density  ratios.  Accordingly,  any  fluctuations  that  are  random,  and  that  may  be 
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present  in  the  ratios  of  the  true  density  to  the  mean  density,  would  tend  to  aver¬ 
age  out  in  computing  the  ballistic  density. 

Another  question  that  deserves  consideration  is— how  important  is  geographic 
location,  season,  and  time  of  day?  Only  five  stations  used  in  the  study  are 
located  in  wide  ocean  areas,  far  from  land,  and  between  30°  and  50*  N  latitude. 

Data  were  available  for  only  two  times  of  day,  0000  and  1200  hr  G.m.t.,  which  makes 
it  difficult  to  determine  the  diurnal  variations.  Based  on  the  data  available, 
there  does  not  seem  to  be  much  payoff  for  preparing  separate  tables  for  different 
times  of  day  and  different  seasons  for  the  geographic  areas  that  might  be  repre¬ 
sented  by  the  five  stations. 

There  does  appear  to  be  a  significant  difference  between  the  correlates  of 
stations  E,  N,  and  C,  North  latitudes  35",  50“ ,  and  52“  45',  respectively.  The 
ballistic  density  at  30,000  ft,  corresponding  to  the  same  surface  density,  is  about 
one-half  percent  higher  at  the  lower  latitude  stations  than  at  station  C.  The  sur¬ 
face  densities  run  higher  at  station  C.  The  upper  air  density  may  decrease  more 
rapidly  than  it  does  at  the  other  stations,  resulting  in  lower  ballistic  densities 
for  the  same  surface  density. 

Data  by  season  do  not  appear  to  correlate  differently,  except  at  station  C. 

The  largest  difference  noted  here  was  for  surface  density  1.03,  where  the  ballistic 
density  at  30,000  ft  was  about  one  percent  lower  in  summer  than  in  other  seasons. 

Wiile  the  procedure  presented  in  this  report  has  been  used  on  only  five  sta¬ 
tions  that  have  an  oceanic  climate,  it  should  be  possible  to  use  essentially  the 
same  procedure  for  processing  data  from  stations  where  the  climate  is  more  affected 
by  land  masses.  Some  modifications  might  be  required,  since  the  temperature  lapse 
rate  might  be  quite  different.  The  diurnal  and  seasonal  variations  in  pressure  and 
temperature  are  normally  much  larger  over  land  masses  than  they  are  over  large 
bodies  of  water. 

A  similar  but  earlier  study^  was  made  by  the  Army  Electronics  Command,  Fort 
Monmouth,  New  Jersey.  Data  from  land-based  stations  were  used  in  this  study  and  a 
number  of  tables  were  prepared  for  each  geographic  area,  season,  and  time  of  day. 
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CONCLUSIONS 

1.  There  does  appear  to  be  a  good  correlation  between  the  surface  density  and 
the  density  at  various  altitudes  over  the  mid-ocean  areas  studied. 

2.  The  resulting  correlations  are  different  for  different  latitudes.  They  do 
not  change  appreciably  with  diurnal  time  or  season  over  the  areas  studied. 

3.  The  differences  between  surface  and  upper  atmospheric  conditions  are  much 
larger  over  littoral  areas. 

4.  The  procedure  described  in  this  report  for  determining  the  correlation 
between  upper  air  density  and  the  surface  density  may  be  used  for  a  variety  of 
stations  with  only  minor  modifications  to  the  procedure. 

5.  While  the  procedure  for  estimating  ballistic  density  was  developed  prima¬ 
rily  for  firing  Navy  guns,  it  may  also  be  useful  to  predict  the  ballistic  den¬ 
sity  for  air-launched  missiles. 

6.  It  is  estimated  that  the  density  ratios  at  altitudes  not  to  exceed 
36,000  ft  may  be  predicted  with  an  error  not  to  exceed  1.0  percent,  using  the 
correlation  between  the  surface  density  and  the  density  of  the  upper  atmo¬ 
sphere. 

7.  It  is  estimated  that  the  ballistic  density  may  be  predicted  with  an  error 
not  to  exceed  0.5  percent  in  90  percent  of  the  cases  where  the  altitude  does 
not  exceed  36,000  ft. 

8.  It  is  estimated  that  the  ballistic  density  may  be  predicted  with  an  error 
between  one  and  two  percent  where  the  altitude  is  between  36,000  and 

50,000  ft.  This  error  could  be  reduced  somewhat  if  the  data  used  in  the  study 
had  been  collected  with  present  day  equipment. 
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RECOMMENDATIONS 

1.  Correlations  should  be  obtained  for  other  latitudes,  including  both  mari¬ 
time  and  beach  areas  of  strategic  importance. 

2.  More  data  should  be  obtained  for  the  extremely  low  and  high  surface  densi¬ 
ties. 

3.  A  further  study  should  be  conducted  to  determine  how  the  procedure  of 
estimating  the  ballistic  density,  as  described  in  this  report,  could  be  imple¬ 
mented  in  fire  control  systems  and  in  preparation  of  range  tables. 

4.  Consideration  should  be  given  to  using  the  analytic  expression  for  density 
in  computer  programs  that  simulate  the  trajectories  of  unguided  missiles. 
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A  NONLINEAR  LEAST-SQUARES  TECHNIQUE 


Statement  of  the  Problem 


Let  (x^,  y.) ,  i  =  1,2,  ...,  N  be  a  set  of  N  data  points  and 

f(x,  3^  |  a,^  >  • » » >  fljj)  C  A— 1 ) 

be  a  function  with  x  as  the  independent  variable,  (aj),  j  =  1,2,  ...,  M  are 
parameters  of  the  function  and  the  function  may  not  be  linear  in  these  parameters. 


We  would  like  to  determine  values  of  the  M  parameters  such  that 


i=N 

F(a^,  a^,  a^)  £  ^i^i  —  a£’  a2* 


(A-2) 


is  a  minimum.  This  is  the  least-squares  solution,  which  is  the  same  as  minimizing 
the  sum  of  the  squares  of  the  residuals.  A  weight  W^  is  assigned  to  each  data 
point . 

For  simplicity  of  notation,  we  let 

f^a)  =  fCx^,  aj,  a2>  ....  ,  i  *  1,2,  ....  N  (A-3) 

where  a  =  (a.,  a0 ,  ...,  au) .  The  only  restriction  on  the  function  f.(a)  is  that 
1  JL  M  1 

the  partial  derivatives  of  the  function  with  respect  to  the  M  parameters  be 
continuous  in  all  of  the  arguments. 

This  implies  that  f^(a)  is  a  function  that  can  be  differentiated. 


We  assume  that  these  partial  derivatives  exist  at  a  point  a,  and  in  some 
neighborhood  about  I,  which  includes  the  point  a.  Let 

a  ■  a  +  Aa  (A-4) 
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where  a  -  (5^ ,  a 2>  aM>  and  Aa  *  (Aa1 ,  Aa2>  ....  Aa^. 

If  f.(a)  can  be  differentiated  at  a,  then  the  differential  at  a  is  an 
1 

approximation  to  f^(a)  -  f^(I)  and  f^(a)  -  f^(a)  =  0(&a)(a^a). 

Thus , 

j=M  3f.(a) 

f.(a)  -  f.(a)  =  V  — -  Aa.  +  E(Aa)  (A-5) 

1  1  j“i  3aj  J 

where  E(Aa)->0  as  Aa-*0  for  i  *  1,2,  ....  N. 


If  the  approximation  a  is  close  enough  to  the  solution  a,  then  the  problem  is 
considered  solved  when 


I  aj  “  5j  M  Aaj  |  <  T.  j  -  1,2,  ....  M  (A-6 ) 

where  T  *  desired  precision  for  the  solution. 

Formulation  of  Solution 


We  now  derive  the  formulation  for  finding  Aa.  The  solution  a  is  then  equal  to 
a  ♦  Aa. 

Using  equation  (A-5),  we  form  the  overdetermined  system  of  equations  (N>M) 


Y. 

l 


f.(a) 


fi<a)  +  l 


j"M  3f.(a) 

— i -  Aa. 

3a.  i 

j-1  j  J 


(A-7  ) 


for  data  points  i  *  1,2,  ...,  N. 
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The  system  of  equations  used  to  find  Aa  becomes 


3f.(a)  3f.(a) 

__i -  Aa  +  — -  Aa0 

3a^  1  2 


3f .(a) 

“5S7"  AaM  =  Yi  -  fi(a) 


(A-8) 


(i  =  1,2,  ....  N). 

The  system  of  observational  equations  (A-8)  can  now  be  written  in  matrix  form 
as  given  by 


where , 


AX  =  Y 


(A-9  ) 


3^(1) 


3fj  (a) 


3^(1) 


3f2(a) 


3f2(D 


«2(i> 


V5> 


V5) 


3V« 
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yi  -  fx(«) 
y2  ~  f2(5) 


yN  -  fN(i) 


The  column  matrix  Y  is  called  the  residual  matrix.  Solving  matrix  equation 
(A-9)  yields  the  normal  equations  in  matrix  form 


T  T 

A  AX  =  AY 


(A-10) 


If  the  residuals  are  to  be  weighted,  we  can  include  the  weight  matrix  W 
where. 


W 


W, 


w, 


N  J 


(A-l 1 ) 


Completing  the  solution  to  the  matrix  equation  (A-10)  and  upon  including  the 
weight  matrix,  we  have 


X  *  (ATWA)_1  ATWY  (A-12) 

X  is  the  least-squares  solution  to  the  system  of  equations  (A-8)  with  the 
weights  included.  Elements  of  the  X  matrix  are  now  tested  in  absolute  value 
against  the  tolerance  T.  If  the  tolerance  test  is  satisfied,  we  are  done  and  the 
least-square  solution  to  equation  (A-2)  is  given  by 

■  5^  +  Aaj,  j  -  1,2,  ...,  M  (A-l 3) 
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Whenever  the  tolerance  test  fails,  equation  (A-13)  gives  the  next  approxima¬ 
tion  and  the  iteration  process  is  repeated  with  matrix  equation  (A-12)  until  the 
tolerance  test  is  satisfied. 

Remarks 


Once  an  acceptable  solution  has  been  found,  it  may  be  desirable  to  know  how 
well  the  function  fits  the  given  data.  Analysis  for  how  well  the  solution  fits  can 
usually  be  made  by  examining  the  residuals  for  each  observation. 

The  variances  for  the  solution  (a)  can  be  found  from  the  variance,  covariance 
matrix. 

If  the  variances  and  residuals  are  large,  a  different  approximating  function 
f^(a)  may  be  tested. 

When  polynomials  are  used  for  testing,  the  residuals  will  drop  off  rapidly  as 
the  degree  of  the  polynomial  increases  from  one.  Further  increase  in  the  degree  of 
the  polynomial  will  show  little  improvement  of  the  residuals. 
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APPENDIX  B 

AN  ALTERNATE  METHOD  FOR  COMPUTING  DENSITY 
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AN  ALTERNATE  METHOD  FOR  COMPUTING  DENSITY 

Since  the  procedure  used  to  fit  the  observed  densities  provided  little  infor¬ 
mation  as  to  why  the  correction  term  for  altitudes  below  3,500  ft  improved  the  fit, 
an  investigation  was  made  using  a  sample  set  of  data  obtained  from  25  sondes.  The 
sondes  were  taken  consecutively  from  the  group  that  had  been  sorted  according  to 
season  after  the  surface  densities  were  arranged  in  numerical  order.  The  median  of 
the  surface  densities  in  the  sample  was  equal  to  the  U.S.  Standard  density  at  mean 
sea  level.  The  sondes  were  released  at  station  E  during  the  spring  season, 

1200  hr  G.m.t. 

Separate  least-squares  fits  were  made  to  the  pressure  and  temperature  in  the 
altitude  range  from  0  to  36,089  ft.  No  attempt  was  made  to  fit  the  temperature  or 
pressure  above  36,089-ft  altitude. 

The  equation  used  to  fit  the  pressure  is  similar  to  equation  (l) 

(p,h) 

P  *  ?1e  +  p3h  (B-l) 


where 


P  is  the  pressure  (mb) 
h  is  the  altitude  (ft) 

Pj ,  Pj,  and  p^  are  constants  determined  in  the  least-squares  fit 

The  observed  temperatures  were  first  converted  to  virtual  temperatures  with 
the  equation 


where 


TV 


1.0 


T _ 

0.3783*RH*V 

P 


(B-2) 


TV  is  the  virtual  temperature  (*K) 

T  is  the  observed  temperature  (’K) 

RH  is  the  observed  relative  humidity  (dimensionless) 
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V  is  the  saturation  vapor  pressure  (mb) 

F  is  the  observed  pressure  (mb) 

Several  least-squares  fits  were  made  with  the  virtual  temperature  in  the  alti¬ 
tude  range  from  0  to  36,089  ft  using  polynomials  that,  in  every  case,  produced 
large  residuals  below  3,000  ft.  A  big  improvement  was  made  by  segmenting  the  data 
into  three  zones,  0  to  3,500,  3,501  to  6,000,  and  6,001  to  36,089  ft.  A  one-degree 
polynomial  least-squares  fit  was  used  in  the  two  lower  zones  and  a  third-degree 
polynomial  was  used  in  the  upper  zone.  The  intersections  of  the  equations  are 
located  at  3,731  and  5,270  ft  which  were  used  to  define  the  boundaries  of  the 
zones. 

The  density  was  computed  with  the  equation 

p  -  P / ( R*T )  (B-3) 


where 


p  is  the  density  (lb/ft^) 

P  is  the  pressure  determined  by  a  least-squares  fit  (mb) 

T  is  the  temperature  determined  by  a  least-squares  fit  (°K) 

R  is  a  constant,  45.981610868 

Computing  the  density  in  this  manner  made  a  small  reduction  in  the  mean  of  the 
residuals  but  no  improvement  in  the  standard  deviation.  While  a  single  function 
fits  the  pressure  from  sea  level  to  36,089  ft,  the  nonlinearity  in  the  temperature 
profile  is  quite  apparent.  The  lapse  rate  in  the  temperature  based  on  the  least- 
squares  fit  from  sea  level  to  3,731  ft  is  2.60'C  per  1,000  ft  and  0.358°C  per 
1,000  ft  from  3,731  to  5,270  ft  using  the  intersections  of  the  equations  to  define 
the  zone  boundaries.  The  average  lapse  rate  based  on  the  least-squares  fit  in  the 
zone  from  5,270  to  36,089  ft  is  1.99’C  per  1,000  ft.  The  latter  value  is  in  good 
agreement  with  the  lapse  rate  of  1.98*C  per  1,000  ft  from  sea  level  to  36,089  ft 
that  was  used  to  define  the  U.S.  Standard  atmosphere. 

Profiles  of  the  temperature  and  pressure  obtained  with  the  least-squares  fit 
are  shown  in  Figures  B-l  and  B-2,  respectively.  A  profile  of  the  density  computed 
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FIGURE  B-2.  STATION  E  1200  hr  G.m.t. 

Pressure  profile  based  on  a  least-squares  fit  with  an  exponential  function 
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with  equation  (7)  is  shown  in  Figure  B-3.  For  comparison,  Figure  B-4  is  a  profile 
of  the  density  computed  with  equations  (1)  and  (3)  where  the  coefficients  were 
obtained  by  making  a  least-squares  fit  to  the  observed  densities.  Very  little  dif¬ 
ference  can  be  seen  when  curves  of  Figures  B-3  and  B-4  are  superimposed. 

From  these  results— (1)  the  small  difference  in  the  residuals  between  the  two 
methods  of  calculating  densities,  and  (2)  indiscernible  difference  in  the  two 
curves  obtained— it  appears  that  there  is  no  advantage  in  computing  density  from 
smoothed  temperature  and  pressure  data  versus  smoothed  density  data.  It  may  be 
noted  that  the  first  method  requires  11  undetermined  constants,  the  latter 
only  4. 

The  method  of  making  separate  least-squares  fits  to  the  temperature  and 
pressure  does  provide  a  much  clearer  explanation  as  to  why  the  correction  term  in 
equation  (3)  improved  the  fit  to  the  density.  The  reason  is  the  large  variations 
in  the  lapse  rate  of  the  temperature  at  the  lower  altitudes. 

The  equations  used  to  fit  the  temperature  are 
Altitude  (ft) 

0  -  3,500  Temperature  =  a^  +  a  h 

3,500  -  6,000  Temperature  =  a^  +  a|h 

6,000  -  36,000  Temperature  =  a^  +  a"h  +8^^  +  a'.jh'* 


where 


temperature  is  the  virtual  temperature  °K 
h  is  the  altitude  in  feet 

The  altitudes  of  intersections  are  3,731  and  5,270  ft. 
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FIGURE  B-3. 

Density  profile  computed  from  the  data  of  Figures  B-l  and  B-2. 
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FIGURE  B-4.  STATION  E  SPRING  1200  hr  G.m.t. 

Density  profile  computed  directly  from  a  least-squares  fit  to  the  density. 
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The  coefficients  obtained  in  the  least-squares  fit  are 


aQ  =  0.29098019920169  E+3 

a:  =  -0.26010921095663  E-2 
a'Q  =  0.28261034232183  E+3 

aj  =  -0.35814511053088  E-3 
aj^  =  0.28474414363984  E+3 

a^  =  -0.42144325442356  E-3 
a'^  =  -0.69303432938534  E-7 
a'^  =  0.85246672752231  E-12 


The  mean  and  standard  deviations  of  the  residuals  divided  by  the  computed 
values  are 


Altitude 

(ft) 

Mean 

Sigma 

N 

0  - 

3,500 

-0.00000002 

0.00603867 

98 

3,500  - 

6,000 

0.00000000 

0.01019354 

76 

6,000  - 

36,000 

-0.00000008 

0.01654041 

357 

0  - 

36,000 

-0.00000362 

0.01433263 

531 

N  is  the  number  of  data  points  in  the  fit. 


The  temperature  profile  is  shown  in  Figure  B-l. 


The  equation  used  to  fit  the  pressure  is 


(c  h) 

P  =  c1  *  e  +  c^h 


where 


P  is  pressure  in  mb 
h  is  the  altitude  in  feet 

The  coefficients  obtained  in  the  least-squares  fit  are 
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cx  -  0.1025894426895  E+4 

c 2  =*  -0.3512226501974  E-4 
c3  =  -0.1469124958633  E-2 

The  mean  and  standard  deviation  of  the  residuals  divided  by  the  computed 
values  are 

Altitude  (ft)  Mean  Sigma  N 

0  -  36,000  -0.00000270  0.01171363  531 

N  is  the  number  of  data  points  in  the  fit. 

The  pressure  profile  is  shown  in  Figure  B-2. 

The  density  profile  shown  in  Figure  B-3  was  computed  with  the  equation 

p  =  P/(R*t) 


where 


p  is  the  density  (lb/ft-*) 

P  is  the  pressure  (mb) 

T  is  the  virtual  temperature  (°K) 

R  is  a  constant,  45.981610868 

The  mean  and  standard  deviation  of  the  residuals  divided  by  the  computed  value 

are 

Altitude  (ft)  Mean  Sigma  N 

0  -  36,000  0.00009937  0.01129087  531 

The  density  profile  shown  in  Figure  B-4  for  altitudes  below  3,500  ft  was  com¬ 
puted  with  the  equation, 
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(c2h)  .  ,  .  si  3500  -  h 

p  =  V  +  c3h  +  (P0  '  Cl}  1  -  3500 


and  for  altitudes  equal  to  or  greater  than  3,500  ft  was  computed  with  the  equation 


(c2h) 

p  =  c  j  e  +  c^h 


which  is  the  same  procedure  described  previously  in  the  report.  The  coefficients 
used  to  create  the  profile  are 

c :  =  0.076640263730  (lb/ft3) 

c 2  =  -0.304780287061  E-4  (ft-1) 
pQ  =  0.076474  (lb/ft3) 

c 3  =  -0.517405589378  E-7  (lb/ft4) 

The  mean  and  standard  deviation  of  the  residuals  divided  by  the  computed  value 

are 


Altitude  (ft)  Mean  Sigma 


0  -  36,000 


0.00048651 


0.01160715 
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APPENDIX  C 

COEFFICIENTS  OBTAINED  BY  FITTING  OBSERVED  DENSITY 
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In  this  appendix,  four  lines  of  data  appear  for  each  group  of  sondes  that  were 
used  in  the  least-squares  fit. 

Lines  1  and  2: 


Item  Identification 


Description 


1 

2 

3 

4 


N 

JE 

AMEAN 

SIGMA 


Number  of  data  points  that  met  toler¬ 
ances 

Number  of  points  used  in  fit 

Mean  of  the  residuals 

Standard  deviation  of  the  residuals 


Line  1  is  for  the  altitude  range  from  0  to  36,089  ft,  and  Line  2  is  for  the  alti¬ 
tude  range  36,089  to  65,000  ft.  The  first  five  items  in  Lines  3  and  4  are  coef¬ 
ficients  ,  C^,  and  for  equation  (1),  which  were  obtained  by  making  a  least- 
squares  fit  in  the  altitude  range  3,500  to  36,089  ft,  and  coefficients  and 
for  equation  (2),  which  were  obtained  by  making  a  least-squares  fit  in  the  altitude 
range  36,089  to  65,000.  The  coefficients  are  listed  in  the  order  of  their  sub¬ 
scripts.  The  next  item  on  Line  4  is  the  median  of  the  surface  densities  of  the 
sondes  used  in  the  fit.  The  last  item  is  a  counter  and  has  no  significance. 

Data  points  refer  to  the  individual  altitudes  tdiere  temperature,  pressure,  and 
humidity  were  recorded.  The  residuals  were  expressed  as  a  ratio  of  the  residual  to 
the  density  obtained  with  the  least-squares  fit  as  stated  previously.  No  data 
points  were  eliminated  in  the  least-squares  fit;  but  tolerances  on  the  residuals 
were  specified,  and  the  number  of  points  that  met  the  tolerances  were  recorded. 

The  tolerance  used  in  the  lower  altitude  region  was  0.025,  and  in  the  higher  region 
it  was  0.150.  A  few  tests  were  made  to  eliminate  bad  data  points  before  the  least- 
squares  fit  was  performed. 
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DENSITY  RATIO  PROFILES  FOR  STATION  E 


Station  E  is  in  the  mid-Atlantic  east  of  Bermuda.  The  uniformity  and  the 
shape  of  the  density  profiles  indicate  that  the  U.S.  Standard  atmosphere  does  not 
represent  the  atmosphere  over  station  E  too  well  for  the  period  sampled.  Almost 
all  the  curves  above  3,500  ft  are  fishhook  shaped.  Below  3,500  ft,  the  density 
ratios  increase  with  altitude.  The  lapse  rate  is  larger  than  standard,  which  is 
referred  to  in  meteorology  as  a  "superadiabat ."  From  3,500  to  35,000  ft,  the  den¬ 
sity  ratio  curves  are  all  concave,  as  seen  from  the  right.  That  is,  density  to 
about  15,000  ft  decreases  with  altitude  more  than  in  the  standard  atmosphere. 

Above  15,000  ft,  the  density  ratios  increase  with  altitude.  Density  ratios  are  all 
greater  than  one  above  30,000  ft,  and  they  grow  larger,  regardless  of  surface  den¬ 
sity  ratios.  Overall,  when  density  is  two  percent  below  standard  at  the  surface, 
it  averages  three  percent  above  standard  at  35,000  ft.  Above  35,000  ft,  the 
curves  are  more  erratic,  which  is  probably  due  to  the  higher  noise  level  in  the 
observat ions . 

Some  seasonal  differences  are  apparent  in  this  region.  For  example,  in  sum¬ 
mer,  both  at  2100  hr  and  0900  hr  local  time,  the  density  ratios  increase  linearly 
from  35,000  to  45,000  ft  to  a  value  of  11  percent  above  standard.  This  seasonal 
behavior  may  be  associated  with  the  seasonal  change  in  the  height  and  temperature 
of  the  tropopause. 

The  diurnal  differences  are  very  small.  When  the  9  a.m.  and  9  p.m.  curves  are 
superimposed,  the  coincidence  of  the  curves  is  remarkable.  This  could  mean  that 
there  are  no  diurnal  changes  in  density;  or  that,  as  the  surface  temperature  and 
density  change  diurnally,  the  density  aloft  changes  in  a  corresponding  manner,  and 
the  correlation  is  preserved.  It  could  also  be  due  to  the  fact  that  9  a.m.  and 
9  p.m.  temperatures  are  nearly  the  same.  At  stations  N  and  V,  the  sonde  times 
would  likely  correspond  to  greater  differences  in  surface  temperature  and  density. 

The  values  for  surface  density  reflect  the  seasonal  changes  in  temperature. 
Surface  density  class  medians  ranged  from  0.95  to  0.98  in  summer;  in  winter,  they 
-anged  from  0.97  to  1.02. 
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APPENDIX  E 


BALLISTIC  DENSITY  PROFILES  FOR  STATION  E 
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A  METHOD  TO  CORRELATE  THE  UPPER  AIR  DENSITY  WITH  SURFACE  DENSIT-- ETC (U) 
JUL  82  L  J  MCANELLY. 
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Local  time  corresponding  to: 


Station  0000  hr  G.m.t.  1200  hr  G.m.t. 
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38 

min 

D 

21 
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02 

hr 

20 
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(To  be  used  with  the  following  tables.) 
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